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DETERMINATION OF THE DAMPING MOMENT IN YAWING FOR 
TAPERED WINGS WITH PARTIAL-SPAM FLAPS 
By Sidney ll. Harmon 
SUMMARY 

A method for determining the, damping moment in yawing 
for tapered wings with partial -span flaps is presented herein. 
Charts are given for untwisted wings with taper ratios of 
0.25, 0-5°> and 1.00, with aspect ratios from 6 to l6, and 
with center-span flaps extending from 25 to 100 percent of 
the wing semispan. The results are also applicable to tip- 
span flaps extending from 0 to 75 percent of the wing semi- 
span. The calculated damping moment in yawing is compared * 
with experimental results for a rectangular wing with a flap 
having a £pan 60 percent of the wing span, ' 

. INTRODUCTION. 

The calculation by Wifeselsberger of the wing damping 
moment for an untwisted elliptical wing in yawing is summa- , 
ri zed and extended in reference 1 to the case of an untwisted 
rectangular wing. Reference 2 presents the results of cal- 
culations/ for a wide range of taper ratio for untwisted wings 
and also for the special angle-^of- attack distribution that 
results from the deflection of partial-span flaps of constant 
chord ratio when the rest of the span is at zero angle of attack 



The results in reference 2 for the yawing derivative due to 
induced drag, however^ contain llnaccuracies because of the 
omission of an important tern in the formula for the yawing 
moment. Also, as noted in reference 2, the results cannot 
be applied .by simple superposition to the case in which the 
lift is contributed simultaneously by partial-span flaps and 
by the plain portions of the wing. This limitation follows 
from the fact that the damping moment varies a.s x the square of 
the angle -of -attack distribution; hence, separate components 
of the lift distribution have interactions that contribute to 
the resultant value of the yawing derivative. 

The present analysis gives the results of calculations 
for the yawing derivative bC n /b (^p) :for untwisted tapered 
wings with partial-span flaps of constant chord ratio at 
various angles of attack. The results are presented for the 
same range of tapbr ratio as is considered in reference 2 and 
for center-span flaps extending from 25 to 100 percent of 'the 
wing semispan. The results preserited for the center-span 

flaps may be applied to tip-span flaps extending from 0 to 

t 

75 percent of the wing semispan. The computations in the 
present paper do not include the part of the yawing moment 
^contributed by the changes in spanwise profile drag* The 
effect of this factor on the 'yawing derivative is discussed 
in reference 2. 
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SYMBOLS 

yawing-moment coefficient due to spanwise induced- 
drag distribution 
angular velocity in yaw, radians per second 
wing span 

wind velocity along plane of symmetry of airplane 
local relative wind velocity at any section % 
coordinate measured along lateral axis of airplane 
circulation around any section 
section lift coefficient 
wing chord at any section . 

induced angle of attack at any section, radians 
coordinate indicating fixed spanwise position 
normal component of velocity 

yawing moment due to spanwise induced-drag distri- 
bution" (C n q s Sb) 
dynamic pressure at plane of symmetry (^PVs 2 ) 
wing area 
density 

parameter defining spanwise position ( 7 - jy cos 0; 

b b \ 

when 9 = 0, 7 ~ ~p when 9 = T = - 

parameter defining fixed spanwise position 
wing chord at plane of symmetry of airplane* 
slope of section lift curve at plane of symmetry of 
airplane, per radian 
• A n coefficients of Fourier series (see reference 5) 



\i parame 
A aspect ratio 



c n r yawing derivative 



. over-all lift coefficient 
^1 ,K 2* and K 5 proportionality constants used to express induced* 

yawing- moment derivative in terms of its 
component parts 
\ taper ratio, that is, ratio of extended tip chord to 

chord at plane of symmetry 
m Q slope of section lift curve, per radian 
B arbitrary constant 

AC T increment of over-all lift coefficient due to deflec- 

tion of flaps 
b f flap span 

0^ section profile-drag coefficient 

Ach increment of section profile-drag coefficient due 

Q Of 

to deflection of flaps 

AC n increment of yawing derivative due to spanwise 
n p 

changes in profile drag 
Subscripts and superscripts: 

w refers to c^ or c do -distribution resulting from 

wing angle of attack 
f .refers to c^- or c d -distribution resulting from 

deflection of flaps 
f c center-span flaps 
f t tip- span flaps 



METHOD AND ANALYSIS 

The method used in the present analysis for calculating 
the stability derivatives is based on the assumptions out- 
lined in reference 2. Because the angular motions con- 
sidered in the present analysis are small (rb/2V 3 < 0,1), the 
influence of the curvature of the wing wake is considered to 
be negligible for practical purposes* Powers of rb/2V g 
higher than unity are also neglected. 

In yawing motion, the following relationships hold: 

where V is the local relative velocity at a section and y 
varies along the span from £ to 




where is the induced angle of attack at any section y 0 

due to the trailing- vortex system. 




where N is the yawing moment due to the spanwise changes 
in induced drag. 
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It is convenient to make the substitution 



y = ^ 00 s 9 

Then, by differentiating equation (1) with respect to the new 
variable 9 , the downwash equation (2) becomes 



a. 



l r 



&{c z c) 



rb 



de~~ + c * c 2v7 sln 9 



d9 



cos 6 - cos 8, 



Equation (Lj.) Is equivalent to the dovmwash equation (8) given 
in reference 2. 

If the Lotz method for determining the span load distri- 
bution (references 3 and I4.) is followed, the circulation 



r = 



c 7 cV c a m 0 V 



2 / A n sin nG 

d(c l c) 



(5) 



Substituting for c^c < and ^ 
in equation (ij.) and integrating gives 



from equation (5) 



c s n 0 



a i 1+b |./ /^n sin 9 2V B Z_ A n 



cos n8 



LI " 1 

The yawing moment N from equation (3) now becomes 



N = -q, 



S 1 l ' * S 7 



c >s 9 | cos G sin 9 c 
(6) 



Integration of equation' (6) and conversion to the 
nondimensional coefficient C n results in 



rb 

2?; 



.V"(2n .+ l)A n A n+1 
1 ' o 



Y_ 2 ^-n 2 " "2 - + ]T (2n + 2)A n A n+2 



(7) 



where 



c 3 m O s 



It will be noted that the even-numbered A-coef f icients 
in equation (7) are directly proportional to the asymmetric 
lift produced in yawing, These coefficients are therefore 
proportional to rb/2V s and their products or powers greater 
than -unity can be. neglected. It follows that, if the even- 
numbered A-coeff icients have the values for unit H)/2V S , 
equation (7) becomes 



6C 



n 



rb 

2Vc 



= c 



n r 



p, 2 TrA 



A,2 » 



" ~f + H (2n + 2)A n A n+2 
1 



(8) 



In equation (8), each of the A-coef f icients is a linear 
function of the angle-of -attack or c^-distribution. For a 
given wing, therefore, similarly numbered A-coef f icients due 
to any number of" cj-di stributions may be superposed, and the 



L 
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sun of these similarly numbered coefficients may be substi- 
tuted in .equation (8) to give the resultant value for C n 

r 

due to the combined lift distributions* * 
The subsequent analysis considers the case of the 
Cj -distribution caused by the combined distributions of the 
angle of attack and the deflection of partial-span flaps. 
The superscripts w and f refer to the c ^-distributions^, 
due to the angle of attack and to the flap deflection, 
respectively. By corabining these c^-distributions , equa- 
tion (8) becomes 
2 



00 



S~ (2n + 1) (A n w + A n f ) (A n+ l W + An*ly 2n(An w + An')' 



I *\2 " 

- i Al I h 1 ,+ 7 (2n + 2r) (a* .+ A n f ) (A n+2 W +. A n+2 f ) 



(9) 



Equation (9) nay be expanded to the form 



+ ^ (2n- + 2 )A n w A^ 



+ )>A n w A„ f ~ A X W Al f + Z (2JV + 2) (A n w A n+ 2 f + A n+ 2 W A n f ) 

1 1 . - 

+ -£ (2n + 1)A/ A n+1 f + £ 2n(A/) 2 - + £ (2n + 2)A n f A n+2 f 

1.1 ■ 1 x , 



(10) 
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The three groups of terms on the right-hand side of 
equation (10) are equal to the yawing derivative due to the 
w-distribution, to the interaction of the w- and 
f -distributions, and to the f -distribution. It should be 
noted that, in* equation (10), the even-numbered A-coef f icients 
have the values appropriate to unit rb/2V g . In the equa- 
tion, each A-coef ficient is directly proportional to the 
over-all lift coefficient that results from the par-, 

ticular c^-distribution to which this A-coef ficient refers. 
It follows that 

c w c wf r f 

n r 0 u n r u n„ p 

n r C^ 2 L w CLw Ac Lf L w L f L f 

where C n ^ w , C n ^ wf , and are obtained from equation 

(10), and the superscript wf denotes the interaction of 
the w- and f -distributions. 

For a given wing, equation (11) may be written 

Cn r = + KaC^ AC Lf + K5 AC Lf 2 (12) 

v/here K]_, K£> and are constants referring to the yaw- 

ing derivative per unit for the angle-of -attack and 

flap-deflection distributions and for the interaction of 
these distributions. 
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RESULTS AND DISCUSSION 

Theoretical Results 

Calculations were made in the present analysis to deter- 

mine the yawing derivative for the lift distribution that 

results from a combined uniform and symmetrical distribution 

of angle of attack and uniform deflection of partial-span 

flaps of constant chord ratio. The computations were made 

by equation (10). The A-coef f icients for the angle-of- 

attack and flap-deflection distributions were obtained by 

the method given in references 2 to ]+. (See equation (9) 

of reference 2.) The computations assumed a value of 5 ^7 

per radian for the slope of the section lift curve. The 

range of the investigation includes three taper ratios: 0.25, 

0.50, and l.OOj three aspect ratios: 6, 10,, and l6j and 

flaps extending from the wing center to 25 to 100 percent 

of the wing semispan* The plan forms of the wings, which 

have rounded tips, are shown in figure 1 of reference 2. ^ 

The results of the calculations are given in figures 1 to 3. 

Figures 1 and 3 are similar to figures 1J and 12, respectively, 

of reference 2. The results presented in these figures in 

-reference 2 are in error because of the omission of the term 
At 2 

_ — which appears in equation (8) of the present paper. 

The variation of the yawing-derivative factor with 
aspect ratio is shown in figure 1 for taper ratios of O.25, 
O.50, and 1.00, for a uniform and symmetrical angle -of -attack 
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distribution, and for AC L ^ =• 0. Figure 1 shows that the 
magnitude of decraases as the aspect ratio decreases 

and as the taper becomes sharper. 

Figure 2 shows the variation with flap span of 
which give3 the yawing derivative due to the interaction 
of the combined distributions . of angle of attack and flap, 
deflection* The actual computations for the taper ratios 
of 0*25, 0»5°* and 1. 00 were made for flap spans of 5° an <* 
100 percent of the wing span. In order to obtain the correct 
fairing for the curves, similar results were computed for the 
elliptical wing with A ~ 10, for flap spans, of 25, 50, 75, 
87«5> an & 10° percent of the wing span. The data of figure 
2 indicate that the magnitude of decreases with aspect 

ratio and increases with flap span. For flap spans that 
are greater than approximately 50 percent of the wing span, 
the magnitude of K2 increases as the wing becomes less 
tapered^ the increase becomes greater as the flap span is 
increased* For flap spans of approximately 50 percent of 
the wing span, varies only slightly with wing taper. 

It is interesting to note in figure 2 that, for f lap-s of 
approximately J>Q percent of the wing span, the curve fox* 
it-, . reverses sign; this condition indicates a destabilizing 
influence for this factor. 

The variation of K5 with flap span for uniform 
deflection of flaps of constant chord ratio and for Ct - 0 
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is shown in figure 3. The actual computations and fairing 
of the curves were obtained in the same manner as those for 
figure 2. The values for show variations with wing 

taper, aspect ratio, and flap span similar to those for K 2 
in figure 2 . 

The total value of the yawing derivative due to induced 
drag for the combined distributions of angle of attack and 
flap deflection is obtained from 'figures 1 to 3 and equa- 
tion (12). 

It should be noted from the derivation of equation (12) 

that the values for C obtained from this equation and 

r 

also from figures 1 to 5 are equally applicable when the angle 
of attack or the flap deflection is either positive or negative, 
.as' long as proper account is taken of the signs* Bec;aiise 
the analysis in deriving the formula for C n shows that , 
K^, and are negative for either a positive or a negative 

angle of attack or flap deflection, the application of 
equation (12) and figures 1 to 3 requires only a considera- 
tion of the sign -of Cj^ r and AC Lf . 

Figures 1 to 5, aa noted previously, are based on a 
value of 5.67 per radian for the slope of the section lift 
curve. ' In order to apply the figures for other values of 
m Q - for example, m 0 = 5»6?B, where B is an arbitrary 
constant - each ordinate is multiplied by B 2 or each 
curve in the figures can be shifted to an equivalent aspect 
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ratio BA, where A is the desired aspect ratio. The 
figures refer specifically to center-span flaps but the 
data may be applied to tip-span flaps of 0 to 75 percent of 
the wing span by the following relationship: 

C nr v+f * = n (c^ + CLfc+ ^ 2 - Kg ( CLtf + c Lfc+t ) A CLfc + K3AC Lfc 2 (13) 

In the equation f^ refers to the tip flaps, f c refers to 

. / ' b f fc \ 
center flaps for which, the' span is (l o 00 - Jb, if 

— is the ratio of the span of the tit) flap to the wing 
b 

span, and f^ +c indicates full-span flaps • The flap deflec- 
tion corresponding to Ct> and Ct> is the same as 
that for the tip flaps* Equation (13), of course, can be 
used in conjunction with figures 1 to 3 only for the 
particular case of the combined ^-distribution that results 
from a uniform deflection of flaps of constant chord ratio 

and from a uniform and symmetrical angle-of -attack distribu- 
tion. . 

Comparison of Theoretical with Wind-Tunnel Results 

The calculated results presented in figures 1 to 3 are 

compared in figure 4 with wind-tunnel results obtained from 

reference 5 for a rectangular wing with part3,al-span flaps of 

constant chcrd ratio. The wing hus the following characteristics: 

Tips Square 

Aspect ratio, A ........... ... 6 

Taper ratio* X ...... L.QO 

bf/b * <%60 

AC Lf ........... 0,56 

c d 0w > P er unit wing span ... 0.024 

Ac<} 0 £, jper unit flap span • 0*08 



14 



and 

C L = C L - AC L 

where C L is the over-all lift coefficient. Prom figures 1 
to 3, based on rounded tips, Ki = -0.0225, K 2 = -0.0219, 
and Kj = -0.0125; therefore, the yawing derivative due to 
induced drag 

c n r = -°- 022 5Cl w 2 - °- 02 19C Lw AC Lf - 0.0125 AC Lf 2 (l[^) 

The yawing clerivative- due to profile drag from equation 
(12) of reference 2 is 

AC np = -0.33(0.02i|) - 0.072(0.08) (15) 

By combining equations (li|) and (15)> the yawing deri- 
vative for the winfe in terms of the over-all ? lift coefficient 
becomes , 

c n r = -0.0225C L 2 + 0.0129l+C L - 0.0177SL 

The comparison of the theoretical values with the wind- 
tunnel results is shown for a range of lift coefficient from 
0 to 1.6. (See fig. !+. ) It will be seen that the agree- 
ment generally is good except at high lift coefficients. 
The greater experimental damping moment at the high lift 
coefficients is probably due to the localized tip vortex 
caused by the square tips, which adds an increment of drag 
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increasing approximately as the square of the lift coeffi- 
cient. 

Langley Memorial Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. 

j 

■»* 
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Figure !•- Variation with aspect ratio of yawing-derivative fac- 
tor due to induced drag for uniform and symmetrical 
angle- of -at tack distribution. ACj, • 0, C n - 2. 
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Figure Zr Variation with flap span of /a*inj-deri votive factor due to induced droa 
for interact/on of uniform and symmetrica/ anak-efaffacK distribution and 
uniform deflection of constant -chord- ratio flaps. c n n = ^C^bC^. 
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Figure J^VarioHon with -flap sport o* yawing-derivotive factor due f 0 
induced draq for uniform def/ecfion of cihstant-chorcl-rotto flops.c, o. 
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Figure 4,- Comparsion of theoretical and wind-tunnel results 
for yawing derivate for partial- span flapped wing. 
A * 6; X = 1.00; b^/b = 0.60; AC L = 0.56; and 0^=0^- 0.56. 
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